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The report shows, for simple cases, under what conditions tail-wheel 
shimmy - sometimes observed in rolling of airplanes - may occur. This 
is done by calculation- of the stability limits decisive for rolling, with 
simplifying assumptions. The model tests are described which were per- 
formed for checking of the calculation results obtained and of the assump- 
tions used. 
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1. INTRODUCTION 


In airplanes manufactured by various firms, oscillation phenomena 
on the tail wheels occurred in rolling and led to failures in some cases . 
One is dealing here with self -excited oscillations of the tail wheel and 
of the tail-wheel fork about the swivel axis of the latter, with the air- 
frame participating in these oscillations to a higher or lesser degree 
according to the ratio of the masses and of the moments of inertia. 

The tail-gear assembly should be constructed in such a manner that the 
swiveling part (that is, tail wheel plus tail-wheel fork) under the influ- 
ence of the restoring forces assumes the mean position without oscillations 

*"Beitrag zur Theorie des Spornradf latterns . " Bericht der Focke- 

Wulf Flugzeugbau G.m.b.H., Bremen, Versuchsabteilung. Technische- Bericht e, 
vol. 7 , no. 2 , 1940, pp. 59 - 70 - 
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when it had been deflected from the rolling direction by any force what- 
ever (aperiodic stable motion) . It is also still permissible that the 
swiveling part, after a deflection, passes several times through the zero 
position before it finally remains in the mean position (strongly damped 
oscillation) . However, the rolling will be considerably disturbed when 
the damping of these oscillations becomes too small, zero, or worse, neg- 
ative. In the last case the oscillations build up from, at first, small 
deflections to larger and larger ones until for some deflection a state 
of equilibrium is again established or failure of a structural member 
occurs. 

The investigations described in the present report had the purpose 
of giving an explanation for the origination of these so-called sh immy 
oscillations and of demonstrating in general the influence of the various 
structural characteristics as far as this is possible in a method not 
applied to definite conditions. The numerical investigation required 
several assumptions and simplifications; in order to check their permis- 
sibility, a series of model tests were performed after the theoretical 
treatment, and were compared with the calculation. 


2. STABILITY CONDITIONS FOR INFINITELY LARGE AIRPLANE MASS 


For clarification of the question under what conditions tail -wheel 
shimmy occurs, the differential equations of motion were set up and 
Hurwitz' criteria of stability were used to examine under what conditions 
these equations have solutions with negative damping, that is, with self- 
excitation. As a result, one obtains one or more equations which indicate 
the boundary between stability and instability. They contain, in general, 
the conditions for static and dynamic stability; if a system is recognized 
as unstable on the basis of Hurwitz' stability criteria, periodic motions 
(tail-wheel shimmy) as well as aperiodic ones may occur. 

In order to minimize nonessentials, the mass of the airplane was 
first assumed to be so large in comparison to the mass of the tail gear 
that the motion of the airplane is not influenced by the oscillations of 
the tail wheel. In the further course of the calculation, it will be 
indicated how the shimmy conditions vary when a motion of the airplane 
in one or several degrees of freedom is considered. 

A number of forces act on the swiveling part which must be in equi- 
librium among themselves; the following symbols are used for the deriva- 
tion of these equilibrium conditions (fig. l) : 

cp shimmy deflection 

displacement of the center point of contact with respect to the 
wheel center plane, cm 
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a distance of the wheel behind the swivel axis, cm 

F 

P load, on tail gear, kg 

Jq, moment of inertia of the swiveling part, referred to the swivel 

axis, cmkgs^ 

c r | the springing coefficient of the wheel for forces acting perpen- 
dicular to the center plane at the point of contact, kg/cm 

c-pcp the springing coefficient of the wheel for moments which twist 
the wheel around the perpendicular through the center point 
of contact, cmkg 

cjpm the springing coefficient of a tail-wheel restoring spring, cmkg 


Cq, the total springing coefficient opposing rotation of the wheel 

about the vertical (cq> = cpcp + c r cp ) , cmkg 

p the resistance coefficient for the rotation about the swivel axis, 

cmkgs 

v rolling velocity, cm/s 

p friction coefficient between ground and wheel 


The following moments act on the swiveling part, referred to the 
swivel axis : 


1 . The inertia resistance 



2. The restoring moment -cpqjCp. In the case of retractable tail 
wheels, there exists usually a spring which restores the tail gear to 
zero position. Even if one imparts to this spring (by catches and similar 
devices), a nonlinear regularity, it is, nevertheless, possible to indi- 
cate for small deflection angles, a constant springing coefficient cpq>. 

A second restoring moment stems from the torsional elasticity of the 
tire; as long as the friction on the ground is sufficient, the tire 
attempts to restore again every rotation about the vertical, since it 
has contact with the runway not in a point but in a surface. The moment 
originating thereby has the value -c c q)Cp, so that, because of 
c<p = c r cp + cpcp, the total restoring moment becomes -c<pcp. 

3. All damping forces and moments which originate by bearing friction, 
material damping, and the like, are, in first approximation, put 
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proportional to the angular velocity of the swiveling part — . They 

dt 

are designated as the natural damping of the oscillation system. The 
sum of all these moments is assumed to be of the magnitude -p 

4. Between runway and tire there originates a frictional force which 
acts at the point of contact at right angles to the wheel-center plane. 
Without discussing here its magnitude in detail, one may state that it 
laterally deforms the wheel. If one designates this deformation by | 
(fig. l), and the springing coefficient of the wheel opposing this defor- 
mation by c r |, the force is c r |£ and its moment about the swivel axis 

ac r |£ . 

The four moments named above must be in equilibrium at any instant. 
Therefore the equation 

" J q> Zf " CcpCP " P ft + aC r^ = ° (1) 

do 

is valid. 

Beside this equilibrium condition for the moments about the swivel 
axis, a second condition may be stated which asserts that the force c r |i 
named in point k is balanced by the component S (lateral force) of the 
frictional force between runway and tire which is at right angles to the 
wheel center plane. The velocity and force relations on the rolling wheel 
can be seen from figure 2. 

The airplane is rolling at the velocity v; this velocity may be 
resolved into a velocity component in the wheel plane (v cos cp) , and into 
a velocity component normal to the wheel plane (v sin cp) . By rolling off 
the wheel would attain a peripheral velocity v cos cp if no slip would 
exist between wheel and runway. Actually, however, the peripheral velo- 
city is smaller ( v ru) since, due to the elastic deformation of the wheel, 

a sliding in the rolling direction occurs . Thereby, a slip of the rolling 
wheel is produced, of the magnitude 

v cos cp - v RU Av 

O — 1 — *— ■ — — — — 

V COS cp y COS cp 

The wheel point of contact has therefore a total velocity v-gp which 

is composed of Av and v sin cp. The frictional force has the value 
pP with a direction opposed to the velocity v R p so that the component 
S perpendicular to the wheel- center plane becomes 



S = pP cos a 


(3) 
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a may be expressed by cp and s 


tan a' = 


Av 


v sin (p t an cp 
If this value is introduced into equation (3), S becomes 


OO 


S 


-UP .... ... 

1 + (s/tan cp)^ 


(5) 


If one plots S against cp according to this equation, the following 
diagram results (fig. 3)* For the wheel which is entirely free of slip 
(s = 0), the normal force increases immediately to its ultimate value 
pP for the deflection zero; for the wheel locked by brakes (s =1), 
this force increases according to a sine law and attains the value pP 
for cp = n/2. The actual wheel lies, according to rigidity and friction, 
between these two lines; it must be noted that the slip s itself varies 
with cp so that the curve S = f(cp) obtains a shape different from the 
form drawn. For the further investigation, the variation in the range 
of small deflection angles is important. The slope of the S/cp curve 
at the point qp = 0 is obtained from 


dS 

dcp 




0 



(6) 


If one replaces the lowest part of the curve by a straight line, which 
can be done with very good approximation, one finds 


S = Pep = p'Pcp 
s 


(7) 


The friction value thus fixed at p' = p/s must be determined from 
experiments for every wheel, every load, and rolling velocity since p 
as well as s are dependent on these parameters . 

If the swiveling part performs a rotation about the swivel axis with 
the angular velocity dep/dt, there originates at the point of contact a 
lateral velocity adep/dt (fig. *0; if* moreover, the wheel is deformed 
by | with the velocity d£/dt, the total velocity normal to the wheel 
plane becomes 

v sin q> + a^+-r~ = v sin cp + acp ' + | ’ 


(8) 
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The velocity of the point of contact varies thereby in magnitude and 
direction. The direction is given by 


tan a = — ( 9 ) 

v sin cp + aq>' + | ' 

and after introduction of s and cp 


tan a 


s 


tan cp + 


arp ' + | 1 
v cos cp 


( 10 ) 


As can be seen immediately from figure 4, the denominator of this expres- 
sion equals tan 3, where 3 is the angle between the velocity of the 
wheel center (vj^j) and the wheel center plane (sideslip or yaw angle) . 

For transition to small angles, one finds for the case where a lateral 
motion of the point of contact is added to the rolling motion according 
to equation ( 7 ) 


S = ^ P3 = p’Pd = p'pfcp + - ^ + i ( 11 ) 

s K V v dt v dt/ ’ 

Since, as was mentioned before, the force S thus found must balance 
the elastic force c r g£, the second condition of equilibrium becomes 

Cr|i + >*'r(<P + \ S + \ S) - 0 (12) 

The two equations (l) and (12) represent the differential equations for 
the motion of the swiveling part and of the point of contact. 

Without discussing in detail the complete solution of the equations 
of motion (l) and (12), one can determine, with the aid of Hurwitz' sta- 
bility criteria, under what conditions the occurring oscillations are 
damped, undamped, or excited. 

If one introduces the expressions 


cp = $e r ^ and | = ~;e r ^ (13) 

wherein r = d + i2n:f (d . . . damping, f . . . frequency), into the 
equations named above and cancels out e 1 ”^, one obtains from (l) 
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(~V 2 " °cp " P r ^ 


+ ac r |H = 0 


m 


= fx;cMi (12) 


(c r | + ^ r)s + n'p(l + * r)® = 0 


(15) 


By elimination of $ and H, there results from (l4) and (15) the main 
equation of the system 

J T i^r 5 + ^Jq)C r g + P Cep + ~y~ a2c rg + P c riJ r + 


t (cgpCrg + ^ , Pac r g) ) 


= 0 


(16) 


It can he shown that the system is stable under the conditions (and only 
under those conditions), that the coefficients Aq to A 3 and the deter- 


minant D = 


a 2 a 5 
A 0 A 1 


are positive. Since the coefficients of the main 


equation always are larger than zero, the stability of the oscillation 
diagram depends solely on the sign of the determinant. If one introduces 
into the latter the values for Aq to A3, one obtains the following 


condition for stable rolling 


p'P< ac r g 


+ P A i 'P c cp + a2 °rg + 
ac r£ ' v *^cp 


Pc r g 

J 9 


+ c. 


r|‘ 


H'P/ 


(17) 


The frequency of the occurring oscillations may be given from the equa- 
tion (2itf )2 = A 2 /Aq as 


ll/cqj + a2c r| + pc r |v/pi'P 

“ 2rt V Jqj 


(18) 


The equation for the oscillation frequency applies rigorously only to the 
undamped oscillation; the self -excited as well as the damped oscillations 
more or less deviate in their frequency from the value thus calculated. 
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From equation ( 17 ) , one recognizes that the rolling velocity v 
appears only in the expression which is dependent on p. If the natural 
damping of the oscillation system is small so that it may be neglected 
with respect to the other forces acting on the tail wheel, the conditional 
equation for stable running assumes the following extraordinarily simple 
form 


p'P<ac r g (19) 

Thus the wheel runs stably as long as the frictional force n’Rp, produced 
by a deflection cp, is smaller than the deformation force axpc r g, occur- 
ring from the same deflection. 

If the complete equation ( 17 ) is solved with respect to p, p indi- 
cates a measure for the magnitude of the required damping if for struc- 
tural reasons the rearward position resulting from equation ( 19 ) cannot 
be adhered to. This case is important particularly for nose wheels. 

One recognizes from equation ( 17 ) that, in the case of swivel axis 
with high friction or with artificial damping, a stable region exists 
even for very small rearward positions since the second term on the right 
side of the equation then becomes larger due to p/ac r g. 

The influence of the rolling velocity will be expounded in more 
detail in the discussion of the test results. 


Tests 

As mentioned at the beginning, the results derived in the previous 
section were checked by means of model tests. For this purpose we used 
the test arrangement described below and represented in figure 5* Since 
it was important to make experiments also for the cases where the air- 
plane mass can no longer be assumed as infinitely large or completely 
rigid, the arrangement of the test apparatus was already provided for 
these conditions; in particular, the rotation of the airplane fuselage 
about the longitudinal axis of the airplane had been taken into 
consideration. 

The pipe 1 represents the mass of the fuselage; it may rotate about 
the longitudinal axis in a self -alining ball bearing 2 and about this 
bearing as point of rotation in the perpendicular (x - z) plane. A motion 
out of this plane is prevented by the bearing-equipped guide 3 . The 
rotation about the x axis is opposed by two springs 4- . Adjustable 
masses 5 are rigidly connected to the pipe by which the moment of inertia 
about the longitudinal axis may be varied within broad limits. The 
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bearing load P of the tail wheel on the runway is adjusted by means of 
the sliding weight 6 and the spring 7- The tail wheel 8 consists of 
soft solid rubber; it is rotatable about the swivel axis 9j its distance 
behind the swivel axis a is adjustable. For the mounting of the swiv- 
elTng part as well as of the tail wheel, ball bearings were used. In a few 
tests, a tail-wheel-restoring spring 10 was installed which restores the 
swiveling part to mean position in the range of small angles with linearly 
increasing moment. The tail wheel runs on a slightly roughened leather 
belt 11 which is driven by an electric motor. 

Since for these first investigations the airplane mass was assumed 

to be so large that it was not influenced by the relatively quick motions 

of the tail gear, the rotation of the pipe 1 about the x axis had to be 
prevented. For this purpose the springs 4 were replaced by rigid con- 
nections and the swivel-axis bearing was laterally restrained. 

The first measurements were to yield the initial values (c r |, c^, 

cptp, Jqn I 1 ') used in the equations (17) and (l8) . The springing coef- 

ficients of the wheel could be determined from the natural frequency of 
the swiveling part. With the belt stationery, the wheel was given a 
slight impact for various wheel-rear positions a and for several tail- 
wheel loads P, and the deflection curve was plotted. If one disregards 
the slight influence of the natural damping p, one obtains 

4rt 2 f 2 = (cq, + a2cr|)/jq,. If one plots therefore against a 2 , 

a straight line must result which has a slope of c r ^ with respect to 
the abscissa and intersects the ordinate scale at the value Cq>. As an 

example, this curve has been plotted for the test series A2 in figure 6; 
since for this series a tail-wheel restoring spring had been installed, 
the oscillation frequency was measured once with this spring and once 
without it, and the springing coefficient cpq, was determined from the 
difference. 

In order to obtain the friction coefficients p'P, the deflection 
forces S were measured normal to the wheel plane for increasing deflec- 
tion angles by means of a spring balance for all tail-wheel loads con- 
cerned; the rearward position of the wheel was selected as far to the 
rear as possible in order to avoid oscillations. The p'P values were 
determined for various rolling velocities; however, the influence of the 
velocity was so small that it remained within the experimental scatter. 

The result of these measurements in the range of small angles (approxi- 
mately up to 5°) is given in tables 3 and 4. 

The moment of inertia of the swiveling part J<p ,increases with 
increasing distance behind the swivel axis a. Its magnitude resulted 
from oscillation tests about the swivel axis set up horizontally for this 
purpose. 
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a = 0.5 1 1.5 2 2.5 3 cm 

Jcp X 10 2 = 4.14 4.27 4.35 4.47 4.59 ^.72 cmkgs2 

4 

Shimmy tests . - In these tests the tail-wheel load P was adjusted, 
for a certain rearward position of the wheel a, in turn to P = 1, 2, 

2.8, and 3*6 kg by shifting of the weight 6 and by compression of the 
spring 7, and then the number rpm of revolutions of the driving shaft 
was regulated by stages between n = 50 and n = 1000/minute corresponding 
to a velocity v = 4l.5 and 830 cm/s. The range of revolution below 
100/minute could not be satisfactorily evaluated because of uneven pro- 
pulsion; the same was true for the range about 600/minute because of 
resonance between the natural frequency of the belt and the propulsion 
rotational frequency. At each rpm-stage the wheel was given a short 
impact; it was then observed whether oscillations of constant or of 
increasing amplitude occurred. If all decisive coefficients (for instance, 
for springing, friction, etc.) were constant - as they had to be assumed 
to be in the calculation - the amplitude of the occurring shimmy oscil- 
lation would increase more and more with time without attaining a fixed 
ultimate value. Actually, however, this is not the case; rather, the 
values are constant only in the range of small angles but vary with 
increasing deflection as was shown, for instance, for the value p'P in 
figure 3- Hence, there appears in many cases a constant deflection which 
is the larger, the smaller the total damping of the system is under the 
various test conditions . 

For a deflection of <J> = ±45°, the swiveling part strikes against 
its limiting stops so that no larger deflections could be measured. 

The tests were subdivided into two series, A1 and A 2, which differ 
by the fact that for the tests of the first series no tail-wheel restoring 
spring was used whereas for the tests of the second series, a spring with 
a restoring moment of 19-7 cmkg had been installed. Due to changes on 
the wheel, the springing coefficients of the latter were, moreover, some- 
what increased. 

The numerical tables 1 and 2 contain the test results of both series 
the deflections and frequencies measured for the various tail-wheel loads,] | 

rearward positions, and numbers of revolutions ^v c m/ s _ 0 . 83n m ^ n "*") . As 
can be seen from the tables, the rearward position was changed each time j|| 
by 0.5 cm; the limiting rearward position for which shimmying stops is, j 

therefore, known, in general, only within this accuracy. When more exact 
data are given below, they have been derived from the variation of the 
deflections as a function of the rearward position. 

.i’ 

The numerical tables 3 and 4 indicate the critical springing and 
friction coefficients for the two series. Furthermore, the result of 
the calculation is given as it was obtained from the equation (19) 
agj. = p'P/c r £. For comparison there appears underneath the limiting 
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rearward position a™ from the testj this is that rearward position of 
the tail wheel for which, in the entire range of revolution (or velocity) 
investigated, no oscillations of constant or increasing amplitude occurred. 

The comparison of test and calculation shows a satisfactory agreement 
so that one may regard as permissible the assumptions made for the inves- 
tigation considered. 


In the model test the rolling velocity did not play an essential role 
for answering the question whether or not shimmy is possible for certain 
tail-wheel loads, rearward positions, springing and friction coefficients. 
If one investigates the terms of the complete equation ( 17 ) 


p m 'P + a^c r ;: pc r | 2 v 

ac ri\ v J cp Jcp + Cr| H'P 

which are dependent on the velocity and on the natural damping, one finds 
that they become infinitely large for v = 0 and v = 00 so that for very 
small and for very large velocities, no shimmy oscillations are possible. 
In between a mini mum value of the expression (20) lies at 




( 21 ) 


For checking of these relations in the model test, we plotted for a 
few tests the variation of the deflections for equal rearward position 
and equal tail-wheel load, but for different velocities, against the rpm 
of the shaft. As an example, we show this for the test of the series A2 
with P = 2.8 kg and a = 2.5 cm (fig. 7)- In this figure the value of 
the expression (20) is plotted neglecting the term Pc r |/<Tcp for the 
values indicated in table 4, with consideration of the ratio between 
velocity and rpm of the driving shaft, p was inserted as 0.0186 cmkgs; 
since for n = 770 /minute the shimmy deflection becomes zero, the left 
side of the equation ( 17 ) is equal to the right one at this point, and 
one may calculate p from it. One sees from the comparison of this curve 
with the variation of 2$ from the test that the deflections 2$ are 
smallest in the range where the damping term is largest, and vice versa. 
Furthermore, one sees that the minimum value of the damping term lies at 
about n = l4o/minute with 0.52 kg so that the error of the simplified 
equation ( 19 ) amounts to only 


0.52 = 0.52 

ac r £ 2.5 X 4.25 


100 percent =4.9 percent 
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One finds a further qualitative confirmation of equation (21) when 
considering tables 1 and 2. On the one hand, the rpm pertaining to the 
maximum deflection increases within the same series (A1 or A2) with 
increasing rearward position a and with increasing tail-wheel load P. 

On the other hand, these velocities lie, for otherwise approximately 
equal conditions, for the series A2 higher than for the series A1 as 
is to be expected on the basis of equation (21) from the difference of 
the c <p values. 

According to equation ( 17 ), it should be possible to find a stable 
range even for very small rearward positions since the second term of 
the right part of the equation then becomes very large. For satisfactory 
installation in the model test (small p) this stable position occurs for 
rearward offsets close to zero, and could therefore not be proved reliably. 


3. INFLUENCE OF THE AIRPLANE MASS 


Whereas, it was assumed in the investigations described so far that 
the airplane, due to its size, is not influenced in its motion by the 
tail-wheel oscillations, it will be indicated in the following section 
how the stability conditions vary if this simplifying assumption is 
dropped. We shall not attempt to give a complete solution for it since - 
with consideration of the numerous degrees of freedom of the motion of 
the airplane - the stability condition would contain so many parameters, 
and in very complicated relations, that a survey of the influence of the 
individual qualifying factors would not be possible after all. Conditions 
become more favorable only where one deals with the complete calculation 
of a definite case for which the majority of structural characteristics 
is fixed and only those need be considered which can still be modified 
by the type of construction of the tail-gear assembly. 

In this part of the report we shall treat, as an example, the sta- 
bility condition for the case that the airplane can rotate elastically 
about its longitudinal axis . Such a motion comes about when the airplane 
oscillates on the two front wheels due to alternating yielding of their 
tires. A similar motion occurs in the case of unsymmetrical forms of 
oscillation of the wing since here also the fuselage (and therewith the 
swivel axis) must perform torsional oscillations about the longitudinal 
axis (x axis) of the airplane. 

For the numerical treatment of the degree of freedom of the motion 
thus fixed, the following symbols are used in addition to those enumera- 
ted in section 2 (fig. 8): 

A deflection angle about the x axis (tramping angle) 

h distance of x axis from runway, cm 
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moment of inertia of the airplane (or of the fuselage) about the 
x axis, 'cmkgs2 

centrifugal moment of swiveling part referred to the swivel axis 
and the x axis, cmkgs2 

moment of inertia of tail wheel, referred to the tail-wheel 
axis, cmkgs2 

angular velocity of the tail wheel 

springing coefficient for the rotation about the x axis, cmkg 
The equilibrium conditions of this oscillation form read: 

1. Moments about the swivel axis (compare equation (l)) 

-Jqp - c^pcp - p ^ + ac r || - J xz - Jo = 0 (22) 


2. Moments about the x axis 


^ ~ PA §7 + hc r^ " Jxz S’ » 0 


3. Equilibrium between frictional force and deformation force (com- 
pare equation (12)) 

c r t! + p 'Pfcp + | ^ ^ + g — ) = 0 (24) 

r 1 ^ V v dt v dt v dt/ 

The following remarks are to be made regarding these equations : 

To equation (22): The first four terms correspond exactly to equa- 

tion (l) . The fifth represents the mass coupling; if the swiveling 
part is deflected by A, there originates at any arbitrary mass point 

d2>, 

of it an inertia resistance of the magnitude -z dm =-£, if z is the 

dt d 

distance of the point from the x axis. If one designates the distance 
from the swivel axis by x, the moment of this inertia resistance about 

a 2.\ 

the swivel axis becomes -xz dm and for the entire swiveling part 


„ „ - d2A _ T d 2 A 
xz dm — - -J X z ~ 
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The sixth term - Joj — indicates the magnitude of the moment which 

dt 

is required for deflecting the tail wheel, rotating with the velocity a > , 

about the x axis with the angular velocity (gyroscopic coupling) . , 

dt 

To equation (23): This equation has a structure exactly corre- 

sponding to that of the previous equation, and expresses the equilibrium 
of the moments about the x axis. 

To equation (24): The equation (ll) states that the force component 

S of the frictional force pP, which is normal to the wheel plane, is 
proportional to the angle of sideslip -9. The angle of sideslip was 
given by the right part of the equation (ll). Due to the motion of the 
swiveling part about the x axis, there originates at the point of contact 

of the wheel, a further velocity component of the amount h — so that 

CLU 

the angle of sideslip becomes for this case: 

■d = t + ~ — + — — + — — (2 6 ) 

v dt v dt v dt 

For the further development, the terms p J<x> and 

Joi are neglected. Whereas the influence of the first two expressions 
dt 

can be estimated from the results of the first part, the influence of the 
gyroscopic coupling cannot be given immediately. One recognizes, however, 
that these terms also make the stability condition dependent on the veloc- 
ity, and that this deviation must take effect for large velocities, since 
Jod vanishes for v = 0. It seemed useful to neglect the expressions 
mentioned above in the calculation, and to clarify the problem of the 
influence of the velocity by model tests carried out simultaneously. 

For the determination of the stability condition, one uses again an 
expression of the form 

cp = $e r ^, % = He 1 "^, 7\ = Ae r ^ ( 27 ) 1 

i 

1 : 

If one introduces these expressions into the equations (22) to (24) and ** 
cancels out e r 4, there results three equations in r, $, H, and K } 
the denominator determinant of which must disappear 


-Jcpr 2 - 

c cp i 

- J xz r2 

9 

ac r | 



"Jxz r ^ 

9 

1 

OJ 

a 

h> 

ca; 

hc r | 

= 0 

( 28 ) 

p'P^l + 


p 'P — r 

V 

9 

c r| + P' p \ r 
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By solution of this determinant one obtains an equation of the 
5th degree in r of the form 

A Q r5 + A-jT^ + AgT? + A^r 2 + A^r + A^ = 0 (29) 


with the coefficients 


A o 

- (Jcp^A _ 

A 1 

= (JcpJ A - 

A 2 

= (Jq) c -^ + 

a 3 

= (Jcp c A + 

a 4 

= (cq> c A + 


a 5 = ( c cp + H 


J xz 2 )ii'P/v 

^xz^) c r| 

<l\ c cp + Jcp c r£h^ + JA Cr ^ a ^ “ 2 JxzCr|ah)n 1 P/v 
J A c cp + J A^’ Pa " J xz^’Ph)c r £ 
c cp c r ^h 2 + c A c r |a 2 )u’P/v 

P a ) c A c r£ 


According to Hurwitz ' investigations, the oscillation form is stable 
when A]_, A^, and, in addition, the following determinants are larger 

than zero after Aq has been made positive.' 


D 


1 



D 2 = 


0 



°3 


a 4 

A 2. 

A 0 

0 


A 5 0 0 

Aj Aj^ A^ 

A 1 A 2 a 3 
0 A 0 A 1 


(30) 


(31) 


(32) 


If even only one of these determinants becomes smaller than zero, 
self -excited oscillations or generally unstable motions of the system 


1 6 
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occur. A survey of the stability conditions thereby obtained is consid- 
erably facilitated by the introduction of the following (dimensionless) 
quantities : 


C <P 

c r |h 1 2 


C A _ |i'P 

c r |h 2 c r| h 


t = 





(33) 


In order to obtain an idea of the influence of the individual param- 
eters on the stability conditions, we performed initially a calculation 
for t, = 0 . £ = 0 signifies that the swiveling part is statically and 

dynamically balanced with respect to the swivel axis, and that, therefore, 
the center of gravity of every section, laid parallel to the x axis, 
falls on the swivel axis. For the system simplified in this manner, the 
stability limits fixed by the determinants D^_, D 2 j ^3 become 


from 

D x > 0: 

x > y(7 - y) 


from 

D 2 > 0: 

px 2 + 7yx > - y 

Ty 2 + (a - 7 2 )y - ay] 

from 

D 5 > 0: 

x[py - (P + 1)7 

^ [7^ + (a. - 7 2 )y - a0 


for x ^ a/p and y > 0 


( 3*0 

(35) 

(36) 


If one represents the limits determined by (34) to ( 36 ) by y = f (x), 
one obtains three limiting curves as plotted in figure 9- The numbers 1, 
2 , 3 assigned to the curves refer to the determinants Dq_, I> 2 , and D 3 . 
Which one of the three curves indicates the actual limit depends on the 
ratio 7 2 /4 : a/p. Figure 9 is drawn for a/p > 7 2 /4. If, however, 

a/p < 7 2 /4, the lowest part of the limiting curve is modified according 
to figure 9a- The curve given by D 2 has not been plotted because it 
lies between 1 and 3 and is never the stability limit. The system 
appears stable when the first and third conditions are satisfied. 


We shall point out in particular two characteristics of the limiting 
curves for the case £ =0: 


1. For x = 0, that is, for 1^=0 or = 00, one obtains y = 7 ; 

if one introduces again the original designations, this yields: 
a/h = p'P/c r |h or ac r | = p'P, thus the condition found in section 2 

(equation ( 19 ) ) for the stability limit for the case of infinitely large 
airplane mass. This same limiting condition is valid with good approxi- 
mation also for the points with small Jq)/J\ as may be recognized from 
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the variation of the curve. One may also expect for this range that the 
dependence on the rolling velocity which is not contained in the equation 
due to the neglect of the natural damping is similar to the one derived 
previously. 

2. The straight line x = a/p is a stability limit. With the orig- 
inal designations, the equation of this straight line is transformed into 

= St or £® = £ 2 ^ that is, the system is stable only as long as the 

C A J cp J A 

natural frequency of the swiveling part about the swivel axis (with the 
springing c^) is larger than the natural frequency of the airplane or 
of the airplane fuselage about the x axis. Closely below this limit 
(see fig. 9 ) the wheel runs evenly even for very small rearward posi- 
tions whereas closely above it no stable run is attainable even for 
very large rearward positions . 


Of particular interest is the question whether the lateral rigidity 
of the wheel c r | has in all cases such a decisive significance for the 

occurrence of tail-wheel shimmy as was found in the investigation of the 
case of infinite airplane mass (see equation (19)). Since the lateral 
rigidity of the wheel c r ^ is in the denominator of the values a, p, 
7 , all three expressions will become zero for the completely rigid wheel 
(c r £ = «>) . Only the ratios (a/ p, p/ 7 ) maintain their former values. 

The limiting curves thereby vary according to figure 10. For 
x < a/p = C cp/ G A the system is always stable, whereas in the region 

x > a/p, the stability range has increased and has been shifted into the 
range of smaller rearward positions. The horizontal asymptote decreases 


from 


P-+ 1 


7 to 7 /p = 


_ |i 'F 


h. 


The foregoing investigations of this section referred to the fully 
balanced swiveling part. If one eliminates this simplification, the 
stability conditions fixed by the determinants D]_, Dg, and D^, turn 

into the following equations; 


from D x > 0: x(l + 7^ - 2£y) > y{y - y) (37) 

from Dg > 0 : (px + a + 7y - 7^x) (x + 7^x - 2 £xy 

+ y 2 - 7y) + (x - £2x2) (pTy _ a (38) 

- py2) > 0 

from Dj > 0 : (x + 7^x - 2£xy - 7y + y2) Ja, 2 + <xyy 

+ x(p 2 y2 - p 2 yy - P7y + 2ap^y 
+ P 7 ^y 2 - a,P 7 ^ - art)] - (x' - £ 2 x 2 ) 

X (P 7 y - a - py 2 )2 > 0 


( 39 ) 
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If one plots, in the same manner as for the mass -balanced swiveling 
part, y as a function of x, one obtains a variation similar to fig- 
ure 9 ’ figure 11. One notices as an essential difference that the rec- 
tangular set of straight lines ( horizontal straight line from x = 0 
to x = a/p and vertical straight line through x = a/p) has become a 
curved line passing through the origin. In the range of small x values 
the branch of the limiting curve going through x = 0 , y = 7 has not 
been changed greatly. 

In the case of the balanced swiveling part, it has been pointed out 
that the natural frequency of the swiveling part c<p/Jcp must be higher 
than that of the fuselage if stable rolling is to be attained without 
damping with wheel rearward positions equal to or smaller t han n'P/c r |. 
For the nonbalanced swiveling part one finds the corresponding condition 
by searching for the point of the limiting curve for which y = 7 . If 
one substitutes y = 7 into the equation (39 ) > one obtains after a few 
transformations 


x 


a 


P1 + £ (£' r ) 


(4o) 


Thus it is valid in first approximation for this case, too, that one must 
have x < a/p, that is, C cp/5cp > c^/j^ if shimmy oscillations are to be 
prevented for relatively small rearward positions. Figure 12 shows the 
difference of the stability limits for a balanced swiveling part (£ = 0 ) 
for the following relationships 


a = 0.0174 
P = 1.2 
7 = 0.115 

£ = 0 , respectively, 2.25 
Tests 

Adjoining the calculation of the stability conditions, several model 
tests with the arrangement described in section 2 were carried out. They 
are to serve on the one hand for control of the calculation and its 
assumptions, and, on the other hand, are to show the influence of the 
rolling velocity which had dropped out of the stability equations due to 
neglect of the natural damping and the gyroscopic coupling. From the 
various test series which were performed for several tail -wheel loads 
and springing characteristics of the tail wheel and of the airplane fuse- 
lage about the x axis, two will be singled out here: One series with 

balanced swiveling part (series B 7 ), and one with non-mass -balanced 
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swiveling part (series Cl) . For both series we plotted, on the basis of 
the stability conditions derived, the boundary y = f(x) for values a, 

J3, , y kept constant and then adjusted successively different rearward 
positions a of the tail wheel for a certain J y. Thereafter all rolling 

velocities in question were executed for every a; every time the wheel 
was given a short impact, and the occurring oscillations were observed. 

The results are shown in the curves (figs. 13 and Ik) by means of the 
type of designation of the test points. Points with equal moment of 
inertia have a common series symbol (for instance B 71 , B 72 , Cll, 

C12, etc . ) . The characteristic values for the two series may be taken 
from tables 5 and 6. 

The introduction of nhe value £ - which had taken place already 

in the calculation - requires a certain explanation. £ was determined 
by £ = Jxzy^cp- With increasingly rearward position of the wheel, J xz 
increases; so does J^,, but not in the same manner. Therefore, it cannot 
be expected that £ remains equal for all rearward positions. And, it 
is not the case, either. Rather, £ fluctuates for the test series Cl 
mentioned according to the following table 7- 

The comparative calculation should be carried out for all these £ 
values. This was omitted, however; we calculated throughout with 
£ = 2.25. We thus achieved numerically an essential simplification with- 
out significantly falsifying the result. The value 2.25 was singled 
out because it is a good mean value for the mean rearward positions which 
are of the greatest interest. On which side the error is made at the 
individual point may be estimated from consideration of figure 12, since 
here the curves for £ = 0 and £ = 2.25 are plotted. 

Regarding the value £ = 0 of the series B7, it should be noted that 
it was realized in the test arrangement as closely as possible. The 
static equilibrium of the swiveling part about the swivel axis was 
achieved completely, the dynamic equilibrium, however, only with a more 
or less satisfactory approximation. Possibly this error could be the rea- 
son for the deviation between test and calculation in the region of the 
vertical stability limit at x = 0.0232. 

In general, a comparison between test and calculation shows good 
agreement. Whereas, in the series B7, the influence of velocity is insig- 
nificant, it is found in series Cl in the sense that in the range of small 
x values, shimmy occurs particularly at small velocities and ceases in 
case of more rapid rolling, as known from the tests with infinitely large 
airplane mass. In contrast, the influence of the rolling velocity in the 
range of larger x values is exactly the opposite. 
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k. CONCLUSIONS AND SUMMARY 


The conditions for a stable run of the tail wheel, that is, for the 
avoidance of shimmy oscillations were numerically investigated, under 
simplifying assumptions; the results obtained were checked by model tests 
The essential assumption was that the center of gravity of the airplane 
moves rectilinearly in the rolling direction. The following degrees of 
freedom were treated: 

1. Rotation of the swiveling part about the vertical swivel axis 
plus lateral deformation of the wheel, and 

2. Rotation of the swiveling part about the swivel axis plus lateral 
deformation of the wheel plus rotation of the airplane, and therewith of 
the swivel axis, about the longitudinal axis of the airplane. 

The complete calculation of the first case yielded as the most impor 
tant result that the stable run of the tail wheel is dependent on a suit- 
able rigidity of the wheel and, of course, also of its support in the 
lateral direction. In the case of practically frictionless mounting of 
the swivel axis, the equation (19) : p'P < ac r ^ is decisive which states 

that the run of the tail wheel is stable as long as the frictional force 
(normal to the wheel plane) originating for a deflection cp is smaller 
than the deformation force acpc r | appearing for the same deflection. 

For a given tail-wheel load, the stability is thus favorably influenced 
by a large distance of the wheel behind the swivel axis a, a large lat- 
eral rigidity of the wheel c r £, and a small ground friction p ' . The 
rolling velocity has no effect on the onset of the shimmy oscillations. 
If, however, the mounting is affected by friction, or if dampers had to 
be applied for attainment of rolling stability, the stability conditions 
are regulated according to equation ( 17 ) which also shows the dependence 
on the rolling velocity. 

For the second case where the rotation of the swivel axis about the 
longitudinal axis of the airplane was taken into consideration, the same 
distance of the wheel behind the swivel axis a = p'P/c r £ is sufficient 
as long as the natural frequency of the swiveling part about the swivel 
axis - given by c^/jq, - is sufficiently large with respect to the 
natural frequency of the airplane or the fuselage about the longitudinal 
axis of the airplane. This can be concluded from the behavior of the 
swiveling part possessing mass balance where for Cqj/Jq; < c^/j^ sh immy 

occurs for all distances of the wheel behind the swivel axis (fig. 9)- 
For the nonbalanced swiveling part, the exact limit is fixed by the 
equation (hO) . If the ratio of the moments of inertia is larger 

than the value determined by this equation, shimmy can be prevented only 
by very large distances of the wheel behind the swivel axis. Only when 
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the two natural frequencies lie sufficiently far apart, smaller distances 
of the wheel behind the swivel axis are again sufficient. In the range 
mentioned, tail-wheel shimmy may occur also for a perfectly rigid wheel 
(fig. 10 ), since the lateral elasticity of the wheel is replaced by the 
torsional elasticity of the fuselage (c?J. 

It appears expedient to determine, for the tail-wheel loads in ques- 
tion, the values |i ' (for instance, on a, concrete runway), c r | and 

therewith a gj. = p'P/cj,^ of the usual tail-wheel types by experiments 
since one must not select a distance behind the swivel axis smaller than 
agr if one wants to enforce a stable run without damping. The two quan- 
tities are the only qualifying parameters which must be known about the 
tail wheel to make the design of a stable assembly possible. 

All considerations and tests were carried out with the swivel axis 
in vertical position. On the basis of the considerations which led to 
the determination of the stability conditions, it. is immediately possible 
to draw conclusions as to the behavior in case of an inclination of the 
swivel axis, the more so, since for the usual tail-gear assemblies, only 
small angles of inclination have to be considered. 

If one introduces into the differential and stability equations the 
normal distance of the point of contact from the swivel axis a (fig. 15 ) 
as the distance of the wheel behind the swivel axis, all equations remain 
valid, and it is clear that, for equal design of the tail gear, a will 
be the larger and the stability condition the more easily satisfied, the 
steeper the inclination of the swivel axis. It is true that the inclina- 
tion of the swivel axis produces a restoring moment which is to be derived 
from the tail-wheel load and which should be included into the calcula- 
tion of the amount of c^ so that thereby the natural frequency about 
the swivel axis is increased. However, in the individual case, it must 
be examined whether it is not possible to produce the moment required 
for this purpose by a tail-wheel restoring spring of suitable dimensions, 
without reducing the distance a of the wheel behind the swivel axis. 

A complete solution of the problem requires the consideration of all 
degrees of freedom of the motion of the airplane. As long as any motion 
of the overall center of gravity is disregarded, the following motions 
are of particular interest: the rotation of the airplane about the ver- 

tical axis, and the lifting motion of the swivel axis. This last men- 
tioned motion originates through the deflection of the swiveling part 
about the inclined swivel axis . 

The numerical treatment of these two additional degrees of freedom 
does not show any basic difficulties. A further differential equation 
would have to be set up for each one of them so that the denominator 
determinant becomes of the 7 th or 9 th degree. However, the solution of 



22 


MCA TM 1380 


the corresponding Hurwitz determinants results in a calculation expendi- 
ture which can be justified only in special cases . 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABLE II (Series A 2 ) 



a = 3 cm 
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TABLE VI 


Series 


P kg 

Cep cmkg 
ex cmkg 
c r f kg/cm 
'F kg 


B7 

Cl 

2.5 

2.5 

5A 

33.6 

2320 

2320 

7.2 

4.2 

10.4 

10.4 

0.0278 

0.0174 

1.2 

1.2 

0.115 

0.115 

0 

2.25 


TABLE VII 


a 

cm 

£ 

0 

2.05 

1 

2.19 

2 

2.24 

3 

2.29 

4 

2.26 

5 

2.18 

6 

2.08 

7 

1.97 

8 

1.86 

9 

1.74 
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Rolling direction 



Figure 4.- Velocity and force relations for a deflected tail wheel for 

additional lateral velocity. 
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Figure 6.- Determination of the springing coefficients. P = 2.8 kg; 
c r? = 4.72 kg/cm; c rcp = 8 cmkg; Cpq, = 19.7 cmkg. 
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Figure 7. - Influence of the rolling velocity. 
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Figure 12.- Comparison between swiveling parts with and without 

mass balance. 







